1  // 


AD  -  A 1 3 1  911 


UNCLASSIFIED 


ON  DIFFERENTIAL  EQUATIONS  OF  NONLOCAL  ELASTICITY  AND 
SOLUTIONS  OF  SCREW  D..(U)  PRINCETON  UNIV  NJ  DEPT  OF 
CIVIL  ENGINEERING  A  C  ERINGEN  AUG  83  RR-83  SM-10 
N00014-76-C-0240  F/G  12/1 


NL 


OTIC  FILE  COPY 


ON  DIFFERENTIAL  EQUATIONS  OF  NONLOCAL  ELASTICITY 

AND 

SOLUTIONS  OF  SCREW  DISLOCATION  AND  SURFACE  WAVES 


A.  Cemal  Eringen 
PRINCETON  UNIVERSITY 


Technical  Report  No.  58 
Civil  Engng.  Res.  Rep.  No.  83-SM-lO 


Research  Sponsored  by  the 
OFFICE  OF  NAVAL  RESEARCH 
under 

Contract  N00014-76-C-0240  Mod  4 
Task  No.  NR  064-410 


s 


August  1983 


Approved  for  public  release: 
distribution  unlimited 


Reproduction  in  whole  or  in  part  is  permitted 
for  any  purpose  of  the  United  States  Government. 


DT1C 

ELKCT# 

AUG  3  0 1983 


t>5  08  26  065 


BtCURlTV  CLm^lC*TlQw  or  Tml  PAQE  (Whpn  Ow«  fnif.ttfj 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS  1 

BEFORE  COMPLETING  FORM 

*■  «toO«T  NUMOCN  PRINCETON  UNI  V  1  80',T  ‘CCIJWON  NO. 

■TECHNICAL  REPORT  #58 

1  RCClPlCN  T*S  CATALOG  NuMItR 

*  T,Tl,rt-“""*>0N  DIFFERENTIAL  EQUATIONS  OF 
NONLOCAL  ELASTICITY  5  SOLUTIONS  OF  SCREW 
DISLOCATION  §  SURFACE  WAVES 

Technical  Report 

*  PERFORMING  ORG  REPORT  NUMBER 

83-SM-1 0 

7  AUTHOR).; 

A.C.  Eringen 

•  CONTRACT  or  grant  NUMBER'*/ 

N0001 4- 76-002 040  Mod.  4 

•  PERFORMING  ORGANIZATION  name  and  ADDRESS 

PRINCETON  UNIVERSITY 

Princeton,  NJ  08544 

10  PROGRAM  ELEMENT  PROJECT  TASK 
AREA  0  PORK  UNIT  NUMBERS 

NR  064-410 

M.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

OFFICE  OF  NAVAL  RESEARCH  (code  471) 
Arlington,  VA  22217 

U  REPORT  DATE 

Aug.  1983 

IS  NUMBER  OF  PAGES 

32 

M  MONITORING  AGENCY  NAME  4  AOORCSSof  diff'ont  Irom  Contrail, nj  Othco) 

IS  SECURITY  Class  tot  Out  toponj 

unclassified 

IS#  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 

H  DISTRIBUTION  ITATUtNT  (o I  Mopott) 


*T  DISTRIBUTION  If  »T|m£nT  lot  obttrort  ontmrod  tn  Black  20.  U  dttloconi  tram  Boport) 


»•  lU^LEHENTARv  N  £  T  E  S 


'I  REY  WORDS  fCaofinu*  on  rr»»rt*  aid*  if  n*r(«ian  and  Idonul*  It  hloek  riianM'1 

nonlocal  elasticity,  waves,  dislocation,  surface  waves 


10  ABSTRACT  (Cauiinm  an  »•*•»•»  ai#a  It  naraaii r»  itftnnry  Ay  Natl  inil>Q 

Integro-partial  differential  equations  of  the  linear  theory  of 
nonlocal  elasticity  are  reduced  to  singular  partial  differential 
equations  for  a  special  class  of  physically  admissible  kernels. 
Solutions  are  obtained  for  the  screw  dislocation  and  surface  wave  >. 
Experimental  observations  and  atomic  lattice  dynamics  appear  to 
support  the  theoretical  results  very  nicely. 


DD  iwn  1473  coition  or  t  nov  «•  it  ooiolctc 


■c cunity  CL*snric>Tioi*or  tnk  nee  551  ««<•» •*> 


1 


ON  DIFFERENTIAL  EQUATIONS  OF  NONLOCAL  ELASTICITY 

AND 

SOLUTIONS  OF  SCREW  DISLOCATION  AND  SURFACE  WAVES 


A.  Cemal  Eringen 

PRINCETON  UNIVERSITY 
Princeton,  NJ  08544 


ABSTRACT 

Integro-partial  differential  equations  of  the  linear  theory  of 
nonlocal  elasticity  are  reduced  to  singular  partial  differential  equations 
for  a  special  class  of  physically  admissible  kernels.  Solutions  are  obtained 
for  the  screw  dislocation  and  surface  waves.  Experimental  observations  and 
atomic  lattice  dynamics  appear  to  support  the  theoretical  results  very  nicely. 


»  X 
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1.  INTRODUCTION 

In  the  theory  of  nonlocal  elasticity1,2  the  stress  at  a  reference 
point  x  is  considered  to  be  a  functional  of  the  strain  field  at  every  point 
'  in  the  body.  For  homogeneous ,  isotropic  bodies,  the  linear  theory  leads 
to  a  set  of  integro-partial  differential  equations  for  the  displacement  field, 
which  are  generally  difficult  to  solve.  For  a  spacial  class  of  kernels,  these 
equations  are  reducible  to  a  set  of  singular  partial  differential  equations 
for  which  the  literature  is  extensive. 


The  selection  of  the  appropriate  class  of  kernels  is  not  ad  hoc 
but  fairly  general,  based  on  mathematical  conditions  of  admissibility  and 
physical  conditions  of  verifiability.  For  example,  the  dispersion  curves 
available  from  lattice  dynamics  and  phonon  dispersion  experiments  provide 
excellent  testing  on  the  success  of  these  kernels.  Ultimately,  these 
kernels  should  be  expressed  in  terms  of  interatomic  force  potentials  or 
correlation  functions.  Presently,  several  solutions  obtained  for  various 

problems  support  the  theory  advanced  here.  For  example,  the  dispersion 

1.3 

curve,  obtained  for  plane  waves  are  in  excellent  agreement  with  those  of  the 
Bom-KSrman  theory  of  lattice  dynamics.  The  dislocation  core  and  cohesive 
(theoretical)  stress  predicted  by  nonlocal  theory^’6  are  close  to  those  known 
in  the  physics  of  solids.  Moreover,  nonlocal  theory  reduces  to  classical 
(local  theory)  in  the  long  wave-length  limit  and  to  atomic  lattice  dynamics 
in  the  short  wave-length  limit.6  These  and  several  other  considerations 
lead  us  to  the  exciting  prospect  that  by  means  of  nonlocal  elasticity,  excel¬ 
lent  approximation  may  be  provided  for  a  large  class  of  physical  phenomena 
with  characteristic  lengths  ranging  from  microscopic  to  macroscopic 
scales.  This  situation  becomes  specially  promising  in  dealing  with  imperfect 
solids,  dislocations  and  fracture,  since  in  these  cases,  the  internal  (atomic) 
state  of  the  body  is  difficult  to  characterize. 

The  solution  of  nonlocal  elasticity  problems  are  however  difficult 
mathematically,  since  little  is  known  on  the  treatment  of  integro-partial 
differential  equations,  especially  for  mixed  boundary  value  problems.  There¬ 
fore,  the  treatment  of  these  problems  by  means  of  singular  differential  equa¬ 
tions  is  promising.  This  is  bom  out,  at  least,  with  the  treatment  of  two 
problems  here,  namely  the  screw  dislocation  and  Rayleigh  surface  waves. 
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Results  for  both  problems  are  gratifying  in  that  they  are  supported  by 
atomic  lattice  dynamics  and  experiments. 

The  formulation  of  nonlocal  elasticity  in  terms  of  singular 
differential  equations  is  reminiscent  of  other  fields  of  physics  and  it 
brings  unification  of  various  methodology.  These  are  indicated  briefly 
in  Section  3. 


2.  NONLOCAL  ELASTICITY 

In  several  previous  papers,  we  developed  a  theory  of  nonlocal 
elasticity,  cf.  [1  -  4],  According  to  this  theory,  the  stress  at  a 
reference  point  x  in  the  body  depends  not  only  on  the  strains  at  x 
but  also  on  strains  at  all_  other  points  of  the  body.  This  observation  is 
in  accordance  with  atomic  theory  of  lattice  dynamics  and  experimental 
observations  on  phonon  dispersion.  In  the  limit  when  the  effects  of  strains 
at  points  other  than  x  is  neglected,  one  obtains  classical  (local)  theory 
of  elasticity. 

For  homogeneous  and  isotropic  elastic  solids,  the  linear  theory  is 
expressed  by  the  set  of  equations 


(2.1) 

(2.2) 


tkt,k 


+  p(f* 


0 


tk«(i) 


r 

Ot(  I  X  * -X  I  »T  )  okJ,(x')  dv(x') 


°k<(i'!  =  x  err(;'>  5k«  ♦  2u  ek!(i'> 


(2.3) 
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l2-4) 


3u£(x'). 

wk  ■ 


where  t.  ,  P  ,  f  and  u  are  respectively,  the  stress  tensor, 
mass  density,  body  force  density  and  the  displacement  vector  at  a  reference 
point  x  in  the  body,  at  time  t.  o^(x')  is  the  macroscopic  (classical) 
stress  tensor  at  x*  which  is  related  to  the  linear  strain  tensor  e^fx') 
at  any  point  x'  in  the  body  at  time  t  ,  with  X  and  y  being  Lame 
constants.  The  only  difference  between  (2.1)  to  (2.4)  and  the  corresponding 
equations  of  classical  elasticity  is  in  the  constitutive  equations  (2.2) 
which  replaces  Hooke's  law  (2.3)  by  (2.2).  The  volume  integral  in  (2.2)  is 
over  the  region  V  occupied  by  the  body. 

Field  equations  of  nonlocal  elasticity  are  obtained  by  combining 
(2.1)  to  (2.4).  We  substitute  Eq.  (2.2)  into  (2.1)  and  use  the  identity 


3a 

9xk 


k£v~ 


(x') 


[a  ok£(x')]  +  a 


to  obtain 


(2.5) 


31 / 


C*(|x'-x|)  ok£(x')  n£  da(x')  +  j  a(|x'-x!)  ok£jki  dv(x') 

+  p(fil  -  UP  *  0 


Here  the  first  integral, over  the  surface  of  the  body,  represents  the  surface 
stresses, (e.g.  surface  tension).  Consequently,  nonlocal  theory  accounts  for 


aua&clcz  physic*,  an  important  asset  not  included  in  classical  theories. 
Substitution  of  Eqs.  (2.3  and  2.4)  into  (2.5)  gives  the  field  equations 


(2.6) 


“(Ij’-jD  0  u;,r  «kl  +  *  ui,k)]nk  da' 

a(|x'-x|)  [(X+u)  u^>Jlk  +  uu^kk]dv'  +  p(f£  -  u£)  =  0 


where  a  prime  (')  denotes  dependence  on  x'  ,  e.g.  u*  =  u(x'). 

Integral  equations  (2.6)  must  be  solved  to  determine  the  displacement 
field  u(x,t)  under  appropriate  boundary  and  initial  conditions.  Boundary 
and  initial  conditions  involving  the  displacement  and  velocity  fields  are 
identical  to  those  of  the  classical  theory.  Boundary  conditions  involving 
tractions  is  based  on  the  stress  tensor  tk£  ,  not  on  ok„  ,  i.e. 


(2.7) 


tk£nk  t(n)£ 


where 


are  the  prescribed  boundary  tractions. 


3.  DETERMINATION  OF  NONLOCAL  MODULUS 

From  the  structure  of  the  constitutive  equations  (2.2),  it  is  clear 

_3 

that  the  nonlocal  modulus  o(|x'-x|)  has  the  dimension  of  (length) 
Therefore,  it  will  depend  on  a  characteristic  length  ratio  a/£  ,  where 
a  is  an  internal  characteristic  length  (e.g.  lattice  parameter,  granular 
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distance)  and  £  is  an  external  characteristic  length  (e.g.  crack  length, 
wave  length).  We  may  express  a  in  a  more  appropriate  form  as 

(3.1)  a  *  a( | x *  —x { ,t)  ,  t  =  eQ  a/£ 

where  eQ  is  a  constant  appropriate  to  each  material. 

The  nonlocal  modulus  has  the  following  interesting  properties. 

(i)  It  acquires  its  maximum  at  x‘  =  x  attenuating  with  |x'-x|  . 

(ii)  When  r ->  0  ,  a  must  revert  tothe  Dirac  delta  measure  so  that 
classical  elasticity  limit  is  included  in  the  limit  of  vanishing 
internal  characteristic  length. 

(3.2)  lim  a(|x'-x|,t)  =  6(|x'-x|) 

T  -*■  0 

We  therefore  expect  that  a  is  a  delta  sequence. 

(iii)  For  small  internal  characteristic  lengths,  i.e.  when  t  +  1, 
nonlocal  theory  should  approximate  atomic  lattice  dynamics. 

In  fact,  by  discretizing  Eq.  (2.2),  it  can  be  shown  that  equations  of  nonlocal 
elasticity  reverts  to  those  of  atomic  lattice  dynamics6. 

(iv)  By  matching  the  dispersion  curves  of  plane  waves  with  those  of 
atomic  lattice  dynamics  (or  experiments),  we  can  determine  a 
for  a  given  material.  Several  different  forms  have  been  found1 
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The  following  are  some  examples  which  have  found  applications: 

(a)  One-Dimensional  Moduli 

(3.3)  a(|X|.x)  =  ^  (1  -  M)  ,  |x|  <  £t 

=0  *  |x|  >  £x 

(3.4)  ot ( f x |  ,t )  =  e"lxi/£T 

(3.5)  ot(  I x |  ,t )  =  — - —  exp(-x2/£2t) 

l/nT 

(b)  Two-Dimensional  Moduli 

(3.6)  a(|x|,T)  =  (2r£2i2)  1  K0(/xvx/£t) 

where  KQ  is  the  modified  Bessel  function 

2  -1  ? 

(3.7)  a(  | x |  ,t )  =  (ttt£  )  exp(-  x*x  /  £,  t) 

(c)  Three-Dimensional  Moduli 

(3.8)  a(|x|,t)  =  - 1  3/p  exp(-x*x/4t)  ,  t  =  £2t/4 

8(n  r/z 

Qt(  | x [  ,t )  *  (4tt£2!2)  (x*x)’^  exp(-/x*x  /  t£) 


(3.9) 
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We  note  that  Eq.  (3.3)  gives  a  perfect  match  of  the  dispersion  curve  of 

one-dimensional  plane  waves  based  on  the  nonlocal  elasticity  and  the  Bom- 

3 

KSrmin  model  of  the  atomic  lattice  dynamics  .  For  two-dimensional  lattices, 
Eq.  (3.6)  provides  an  excellent  match  with  atomic  dispersion  curves,  with  a 
maximum  error  less  than  1.2%,  Ari  and  Eringen7.  All  other  nonlocal  moduli 
also  provide  excellent  approximation  to  the  atomic  dispersion  curves,  for  a 

O 

choice  of  eg  ,  Eringen  . 


(v)  We  observe  that  all  nonlocal  moduli  given  above  ar  ialized  so 

that  their  integrals  over  the  domain  of  integration  (line,  surface,  volume) 
give  unity.  Moreover,  they  are  all  6-sequence,  i.e.  when  t  0  we 

obtain  the  Dirac  delta  function.  Because  of  this  property, 

nonlocal  elasticity  in  the  limit  t-*0  reverts  to  classical  elasticity  as 
can  be  seen  by  letting  t  -*  0  in  (2.2),  to  obtain  Hookes  law  of  classical 
elasticity. 

We  now  exploit  this  observation  further  by  assuming  that; 

(vi)  a  is  Green's  function  of  a  linear  differential  operator: 

(3.10)  La(|x'-xj,T)  *  6(lx'-x|) 

If  such  an  operator  can  be  found,  then  applying  L  to  (2.2),  we  obtain 


(3.11) 


L  *k£  *  °k£ 


In  particular,  if  L  is  a  differential  operator  with  constant  coefficients, 
then 


(3.12) 


(L  tki),k 


L  t, 


£,k 


and  (2.1)  gives 


(3.13)  ok£jk  +  L(pf£  -  Pu£)  =  0 

In  this  case,  we  have  partial  differential  equations  to  solve,  instead 
of  integro-partial  differential  equations.®  This,  of  course,  provides  a 
great  deal  of  simplicity  over  the  original  equations  (2.6).  In  particular, 
for  static  problems  with  vanishing  body  forces  (or  more  generally,  when 
L  (pf  -pU  )  =  0  ,  we  have  the  classical  equations  of  equilibrium 

(3.14)  °k;,k  =  0 

which  upon  using  (2.3)  gives  Naviers'  equations. 

The  nonlocal  modulus  (3.8)  is  a  Green's  function  which  satisfies  the 
differential  equation 

(3.15)  V2a-  —  =  0  ,  a(x,0)  =  6(x) 


The  fact  that  this  is  the  case  is  well-known  since  the  solution  t(x,t)  of  the 
diffusion  equation 


n 


(3.16) 


,  at 

H  -  at  =  ° 


which,  fur  t=  0  coincide  with  a  given  continuous  function  o(x),  is  given  by 


(3.17) 


t  = 


a( | x 1 -x 1  ,t)  o(x')  dv(x') 
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when  t'  extends  to  infinity  in  all  directions,  cf.  Courant 
Eq.  (3.17)  is  valid  even  when  V  is  finite  and  the  reference  point  x  is 
not  too  close  to  the  boundaries,  since  a  attenuates  rapidly  to  zero  with 
ix‘-x!  . 

Of  course  the  differential  operator  L  may  be  different  than  the 
diffusion  operator.  For  example,  for  Eq.  (3.6),  one  can  see  that 


(3.18) 


i .e. ,  we  have 


L  -  1  -  A272 


(3.19) 


/ i  2, 2_2» . 

( 1  -  t  s.  .  )  t  =  o 


In  fact,  this  result  can  be  justified  by  an  approximation  of  the  atomic 

a  * 

dispersion  relations.  To  see  this,  consider  the  Bom-Karman  model  of 
lattice  dynamics  and  equate  the  expression  of  the  frequency  given  there 
to  that  of  nonlocal  theory  for  plane  waves. 


(3.20) 


a(k)  =  (cu/cuq)2  =  (2/k2a2)(l  -  cos  ka) 


L 
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We  rev/rite  this  expression  in  the  form 


(3.21) 


2(1  -  cos  ka)/k2a2 


1  - 


k2a2 

~vr 


kV 

360 


(1  +  e2a2k2)-1 


p  O 

for  k  a  <<  1  .  This  is  permissible  since  according  to 

Tauberian  theorems,  the  behavior  of  a  function  near  infinity  is  reflected 
in  the  behavior  of  its  Fourier  transform  near  the  origin.  Thus,  this 
approximation  is  tentamount  to  approximating  a(|x|)  for  large  |x]  . 

In  fact,  by  this  process,  the  boundary  of  the  Brillouin  zone  is  thrown  to 
infinity  so  that  a  does  not  have  a  finite  support.  By  matching  a  = 

(1  +  ega2k2)-1  obtained  this  way  with  Eq.  (3=20)  at  the  end  of  the 
Brillouin  zone,  we  obtain  a  curve  which  approximates  the  atomic  dispersion 
curve  quite  well . 

In  Figure  1,  the  comparison  is  made  for  this  matching,  i.e. 

(3.22)  wa/CQ  =  k  a  ( 1  +  egk2a2)”^  (Nonlocal) 

wa/c0  =  2  s in ( ka/2 )  (Lattice  Dynamics) 

The  matching  is  perfect  at  ka = w,  which  leads  to 

(3.23)  eQ  =  0.39 

The  maximum  error  is  of  the  order  of  6 %  in  [  ka [  <  tt  .  Note,  however, 

that  the  group  velocity  at  ka  =  tt  is  badly  off  and  this  is  the  price  we 
have  to  pay* 
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If  we  accept  this  approximation,  Eq.  (3.21)  upon  inversion,  gives 

(3.24)  (1  -  TVv2)a  =  6 ( |  x | )  ;  T*  =  eQ  a 

The  application  of  this  operator  to  Eq.  (3.17)  then,  leads  to  (3.19).  Of 
course,  other  types  of  approximations  are  possible. 

It  is  not  difficult  to  determine  L  for  other  moduli.  In  fact,  for 
Gaussian  kernels,  (3.5)  and  (3.7)  are  similar  to  (3.15). 

(vii)  The  above  considerations  further  suggest  that  Eq.  (3.17)  may  be 
considered  as  the  probablistic  average  of  o  ,  if  a  is  considered  to 
be  a  probability  density  function.  In  fact,  (3.8)  is  non-other  than  the 
Gaussian  density  function.  Such  a  consideration  resembles  the  method  of 
analysis  of  quantum  mechanics  with  the  probability  density  function  satis¬ 
fying  a  diffusion  equation. 

(viii)  A  somewhat  different  interpretation  of  (3.16)  is  made  by 
considering  a  as  a  correlation  function, in  which  case,  it  should  be 
possible  to  determine  a  from  statistical  mechanical  considerations.  Such 
a  study  is  now  underway. 
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4.  SCREW  DISLOCATION  REVISITED 

Consider  a  screw  dislocation,  in  the  sense  of  Volterra,  located 
at  the  plane  x3=0  of  rectangular  coordinates  xk  ,  Fig.  2.  The  displace¬ 
ment  field  has  only  single  component  u-jU^^.t)  .  The  stress  field  is 
determined  by  solving  Eq.  (3.16)  whose  Laplace  transform  with  respect  to  t  is 


(4.1) 


'2  {k*  -  s  fk<:  '  °k« 


where  ok;  has  only  two  now-vanishing  components 


(4.2) 


31 


=  v 


3u. 

Jx 


1 


°32  =  M 


3u. 


In  (4.1),  a  superposed  bar  indicates  the  Laplace  transform  and 
transform  variable. 

The  divergence  of  Eq.  (4.1), upon  using  (2.1)  and  (4.2), 


(4.3) 


>k£,k  "  '5)esi 


For  the  static  case,  Eq.  (4.3)  gives 
(4.4)  V2  u3  *  0 


s  is  the 

gives 


whose  solution  is 

(4.5)  u3  -  £  tan"1(x2/x1) 


where  b  is  the  Burger's  vector.  In  polar  coordinates  (r,6,z)  this  is 
equivalent  to 


15 


(4.6)  u3  *  EiT  6 

(4*7)  °31  “  *^sine  •  °32  £tFcos6 

Carrying  (4.7)  into  (4.1 ),  we  obtain  differential  equations  for  t31 
and  t^  *  whose  general  solutions,  having  proper  syrmetry  with  respect  to  +. e, 
are  given  by 


t31  =  t^p)  sin  e,  t32  *  T2(p)  cos  e  , 


Ta(p)  =  Aa  Kl<p)  +  Ba  Vp)  +  f  •  0*1,2  • 

p  =  Ss  r  ,  C  =  yb/2iTvT 

Here,  I^p)  and  K-j(p)  are  modified  Bessels'  functions  and  A&  and  Bq  are 
arbitrary  constants.  The  stress  field  must  vanish  as  p-*-®  .  This  implies 
that  Ba  =  0.  In  cyclindrical  coordinates  (r,0,z),  the  stress  field  is 
given  by 

t  =  (A^  +  A2)  K^(p)  sine  cose 

(4.11)  _  27  i  * 

tzp  s  (-A1  «in  6  +  A2  cos  6  )  K^p)  +  C  p 

If  we  imagine  the  edge  line  p«0  of  the  dislocation  as  a  limit  e-*-0  of  a 
small  cylindrical  surface  with  radius  r*e  ,  then  t2r  must  vanish  as  e -» 0. 
This,  through  (4.11),  gives  A2*-A^  so  that 


(4.8) 

where 

(4.9) 

(4.10) 


(4.12) 
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V  ■  0  •  J26  ■  -  A1  *l(p>  +  c  B'1 

The  hoop  stress  t20  will  be  regular  at  p«0,  If  and  only  If,  Aj  *  C. 
Consequently, 

<4-’3>  lz*  ■  &  [  I  -  jz  K1<^>3 

<4-'4)  lze  *  #TF  <>  -  e'rJ/4t)  • 

p 

Remembering  t  *  1 S/4  and  in  terms  of  the  definition  of  a  given  in  Ref.  [4], 
2  2 

t*a  /4k  ,  this  result  is  in  agreement  with  that  found  earlier  in  a 
different  way,4  We  observe  that  the  displacement  field  (4.5)  of  a  screw 
dislocation  in  a  nonlocal  elastic  medium  is  the  same  as  in  a  classical 
elasticity,  even  though  stress  field  is  different.  This  is  the  result 
of  the  particular  choice  of  the  kernel.  The  effects  of  anholonomicity  of 
the  dispersion  curve  near  the  boundaries  of  the  Brillouin  zone  and  the 
nonlinear  force  law  are  ignored. 

Had  we  employed  the  operator  (3.18)  instead  of  (3.15),  we  would 
have  obtained 

(4.15)  t20  *  [1  -  K1  (t*/t£)3  ,  all  other  tfc£  ■  0  . 

This  is  also  regular  for  all  r  in  the  interval  0  <_  r  £*>. 

In  non-dimensional  form,  (4.14)  may  be  written  as 

2 

(4  *>)  Te(c)  -  (1  -  e’p  ) 

wh 

2 

(4.17)  T0  «  2ra  t^/wbk  -  P  (1  -  e‘P  )  ,  P  =  kr/a 
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T.(p)  possesses  a  maximum  at  pc  which  is  the  root  of 

2 

(4.18)  1  ♦  2  p2r  •  e  c  ,  p  =  1.1209 

L  w 


and  ‘emax  is  9iven 


(4.19) 


2  P„ 


emax 


i  +  2  p: 


The  plot  of  Te(p)  versus  p  is  shown  in  Fig.  3.  For  a  perfectly  brittle 

crystal,  the  theoretical  shear  strength  t  is  reached  when  t  Q = t  .  Con- 

y  20  v 

sequently. 


(4.20)  ty  =  0.63S2 

g 

In  a  previous  paper,  we  have  shown  that  k  *  1.65  gives  a  perfect  match  in 
the  entire  Brillouin  zone  of  Born-K£rm£n  lattice  with  an  error  less  than 
0.2*.  Using  this  value  and  b/a  *  1/^  for  fee  materials  in  (4.20)  we  find 


(4.21)  ty/w  =  0.12 

This  result  compares  well  with  the  value  0.11  for  A£  (f.c.c.);  W,  a-Fe 

(b.c.c.)  and  0.12  for  Na  Cl,  MgO.  (cf.  Lawn  and  Wilshaw,  p.  160). 

In  examining  Fig.  2  closely,  we  note  that  when  t  * t  a  brittle 

y  zemax 

perfect  crystal  will  rupture.  However,  if  the  crystal  is  ductile  at  this 

point,  dislocations  will  be  produced.  Thus  the  region  around  the  crack  tip 

0  <  P  <  Pc  is  a  dislocation-free  zone,  i.e.  dislocations  will  emerge  at  p=cc 

and  will  pile  up  in  a  region  p _>  p£  .  This  prediction  of  the  present  theory 

11 

is  supported  by  the  recent  observations  made  in  electron  microscopy  . 
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According  to  the  present  theory,  the  rupture  or  dislocation  initia¬ 
tion  then  begins  not.  at  the  center  of  the  dislocation  p«  0  but  at  a  critical 
distance  p*pc>  0,  This  is  against  our  previous  understanding  of  the  dis¬ 
location  and  fracture  mechanism  but  supported  by  experiments. 

Finally,  we  may  wish  to  enquire  whether  these  predictions  are 
greatly  affected  by  the  kernel  chosen.  In  the  case  of  the  kernel  (3.6), 
leading  to  Eq.  (4.15),  the  TJpJ-curve  is  also  shown  in  Fig.  2  with  p„  , 

v  C 

Temax  and  V  9iven 


If  we  write  h = eQa/0. 3995,  Eq.  (4.23)  agrees  with  the  estimate  of  Frenkel 
based  on  an  atomic  model  (cf.  Kelly^,  p.  12).  In  the  case  of  the  kernel 
(3.6),  the  match  of  the  dispersion  curves  were  provided  for  the  value  of 
eQ  =  0.39  with  an  error  less  than  6%  at  a  point  in  the  Brillouin  zone  (Eq. 
3.23).  Using  this  value,  we  again  obtain  t^/y  =  0.12.  Thus,  in  spite  of 
the  difference  in  the  two  curves  in  Fig.  2,  the  resulting  theoretical 
strengths  are  not  too  far  off  from  each  other.  Moreover,  the  location  of 
the  maxima  are  nearly  coincident  so  that  the  origin  of  the  dislocation 
generation  or  rupture  are  predicted  to  be  the  same. 
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5.  RAYLEIGH  SURFACE  WAVES 

1 3 

In  a  previous  paper  ,  we  have  determined  the  phase  velocity  of 
Rayleigh  surface  waves  and  found  that  they  are  dispersive.  In  these  calcu¬ 
lations,  one-dimensional  kernel  (3.3)  was  used  with  nonlocal  effects  taken 
only  in  the  direction  of  boundary  line,  x2  ~  0  in  two-dimensional  medium 

•  Here»  take  advantage  of  the  two-dimensional 
kernel  (3.6)  which  reduces  constitutive  equations  to  the  form  (3.19).  Upon 
taking  the  divergence  of  (3.19)  and  using  (2.1),  (2.3)  and  (2.4)  with  f=0, 
we  obtain 


(5.D  (x*„j  VkJ*  „Vkk-  (i  -AVi;,  -  o 

For  the  plane-strain,  we  introduce  Lamg  potentials  0(x,t)  and  i|i(x,t), 
(x  *  {x1,x2})  by 


(5.2) 


*  !£_+ 

9x«|  8X2 


leading  to 


u2 


9X£  9x*j 


(5.3) 


c2  v2c  -  (i  -  tW)  ; .  o  . 
C2  V2v  -  (1  -  tVv2)  c  =  o 


where  c1  and  c2  are  respectively,  classical  phase  velocities  defined  by 


(5.4) 


C1 


(A  1 

p  ' 
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We  now  try  solutions  of  (5.3)  in  the  forms  of  surface  waves  characterized  by 
t  =  A  exp[-  kv^2  +  i  k  (x-|  -  ct)]  , 

(5.5) 

v  =  B  exp[-  kv2*2  +  ’  k (xg  -  ct)] 

Equations  (5.3)  are  satisfied  if  are  given  by 

(5.6)  v2  =  1  -  (c/c  )2  [1  -  k2T2£2(c/c J2]"1  ,  a=l,2  . 

Using  (5.5)  in  (5.2),  we  obtain  the  displacement  field  and  carrying  (5.2)  into 
the  constitutive  equations  (2.2)  to  (2.5),  we  arrive  at  the  stress  field. 

Thus,  for  example,  at  x2 =  0  we  have 

t22  =  t(c1/c2)2(v12-l)+2]M1A  +  2iv2M2B 

(5.7)  t21  =  -2iv1M1A  +  (1  +v2)M2B 


where 
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M 


1 


(5.8) 


,co 


dx}  K0  {[(xj  -  x^2  +  x^2]1/2/tO 


0  ’  -OO 

•  exp  [-k  v^x^  +  2k(x^  -  ct)] 


#0° 

dx]  KQ{[(x]  -  x^2  +  x^2]1/2/Te} 

•  »oo 


•  exp  [-kv,x‘  +  ik(x^  -  ct)] 


But  t22  and  t21  must  vanish  at  x2 =  0  .  Hence,  we  must  have 
(5.9)  [lc1/c2)2  ('.  1 2-l  )+2] (1  +  ,22)  -  4  vV,2  =  0 


provided  1  0.  This  is  the  Rayleigh  determinant  for  nonlocal  elastic 

surface  waves.  Upon  carrying  (5.6)  into  (5.9),  eventually  we  rearrange  (5.9) 
into  the  form 


(5.10) 


v[ai 


+  a. 


+  a. 


a4 )  = 


»  -  1  „  m+1  2  .  l+4m-3nr  4  .  ,,  ,6 

al  T6  T  c  +  - S -  c  *  m(l-m)c 


I 


where 
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,,  . .  _  1  ,  2m^-m-3  .2  /, „_2,  4 

(5.11)  +  — 2 e  '  )£  * 

a3  =  |  -  m  +  (2-m-m2)c2  , 

=  1-m  , 

2  l-2v 

Y  =  (c/C 2r  ,  C  =  kx£  ,  rn  »  2(1-vT 

and  v  is  Poisson's  ratio. 

It  is  not  difficult  to  verify  that  for  e  =  0  ,  (5.10)  reduces  to  the 

the  classical  Rayleigh  function  whose  roots  are  recorded  in  Ref.  [14]  for 

various  v.  It  is  also  clear  that  the  roots  of  (5.10)  is  a  function  of  e  , 

consequently,  the  Rayleigh  wave  velocity  is  dispersive.  In  Table  1  below, 

we  give  values  of  /y  =  cR/c2  as  a  function  of  c = eQka  ,  for  various 

Poisson's  ratios  v  .  Phase  velocities  cR/c2  versus  c =  eQka  is  plotted 

in  Fig.  4  for  various  values  of  v  .  The  non-dimensional  frequency  ~e/c?  , 

Table  2  and 

versus  e  are  displayed  in/Fig.  5.  To  provide  a  comparison  of  these  re¬ 
sults  with  the  lattice  dynamic  calculations,  we  divide  the  abscissa  and 
ordinate  of  Fig.  4  by  eQ  =  0.31.  Fig.  6  shows  ka/c2  versus  ka  of  nonlocal 

results.  Comparison  with  atomic  lattice  dynamic  calculations  carried  out  by 
15 

Wallis  and  Gazis  for  KCi  is  excellent.  We  have  used  v  =  0.3  to  approxi¬ 
mate  the  KC*-  molecules  with  an  isotropic  solid.  The  lattice  dynamic  cal¬ 
culations  were  made  for  waves  propagating  in  the  (100)  direction  on  a  (001) 
surface. 


L  i 
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TABLE  I 


Phase  Velocity  of  Rayleigh  Surface  Waves 

CR/C2 


e  =  e0ka 

c 

If 

p 

o 

v  =  0. 1 

v  =  0.2 

v  =  0. 3 

v  0.4 

v  =  0. 5 

0 

0.87402 

0.89311 

0.91099 

0.9274 

0.9422 

0.95532 

0.1 

0.86382 

0.88418 

0.90337 

0.92096 

0.9366 

0.94998 

0.2 

0.8356 

0.85902 

0.88162 

0.90242 

0.92042 

0.93446 

0.3 

0.79476 

0.82176 

0.84857 

0.87377 

0.89519 

0.91025 

0.4 

0.74719 

0.77715 

0.80798 

0.83775 

0.86318 

0.8794 

0.5 

0.69757 

0.72944 

0.76336 

0.79724 

0.82664 

0.84415 

0.6 

0.64889 

0.68173 

0.71768 

0.75479 

0.78772 

0.80646 

0.7 

0.60288 

0.63598 

0.67295 

0.71226 

0.74819 

0.76809 

0.8 

0.56044 

0.59312 

0.63039 

0.67105 

0.7093 

0.73021 

0.9 

0.52178 

0.55374 

0.5907 

0.63193 

0.67187 

0.69363 

1.0 

0.4869 

0.51786 

0.55407 

0.59537 

0.63637 

0.65882 

TABLE  II 


Dispersion  of  Rayleigh  Surface  Waves 
wa/c^ 


e  =  eQka 

< 

11 

o 

• 

o 

v«  0.1 

v  «=  0.2 

v  E  0. 3 

o 

ii 

v  =  0.5 

0 

0 

0 

0 

0 

0 

0 

0.1 

0.086382 

0.08842 

0.09034 

0.09210 

0.09366 

0.09500 

0.2 

0.16712 

0.1718 

0.17632 

0.18048 

0.18408 

0.18689 

0.3 

0.23843 

0.24653 

0.25457 

0.26213 

0.26856 

0.27307 

0.4 

0.29888 

0.31086 

0.32319 

0.3351 

0.34527 

0.35176 

0.5 

0.34878 

0.36472 

0.38168 

0.39862 

0.41332 

0.42207 

0.6 

0.38933 

0.40904 

0.43061 

0.45288 

0.47263 

0.48388 

0.7 

0.42202 

0.44518 

0.47106 

0.49858 

0.52373 

0.53766 

0.8 

0.44835 

0.4745 

0.50431 

0.53684 

0.56744 

0.58417 

0.9 

0.4698 

0.49837 

0.53163 

0.56874 

0.60469 

0.62427 

1.0 

0.4869 

0.51786 

0.55407 

0.59537 

0.63637 

0.65882 
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Air  Force  (Con't.) 

Chief,  Civil  Engineering  Branch 
WLRC,  Research  Division 
Air  Force  Weapons  Laboratory 
Kirtland  Air  Force  Bese 
Albuquerque,  New  Mexico  87117 

Air  Force  Office  of  Scientific  Research 
Bolling  Air  Force  Base 
Washington,  D.C.  20332 
Attn:  Mechanics  Division 

Department  of  the  Air  Force 
Air  University  Library 
Maxwell  Air  Force  Base 
Montgomery,  Alabama  36112 

Other  Government  Activities 

Commandant 

Chief,  Testing  and  Development  Division 
U.S.  Coast  Guard 
1300  E  Street,  NW. 

Washington,  D.C.  20226 

Technical  Director 
Marine  Corps  Development 
and  Education  Command 
Quantico,  Virginia  22134 

Director  Defense  Research 
and  Engineering 
Technical  Library 
Room  3C128 
The  Pentagon 
Washington,  D.C.  20301 

Dr.  M.  Caus 

National  Science  Foundation 
Environmental  Research  Division 
Washington,  D.C.  20550 

Library  of  Congress 

Science  and  Technology  Division 

Washington,  D.C.  20540 

Director 

Defense  Nuclear  Agency 
Washington,  D.C.  20305 
Attn:  SPSS 
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Other  Government  Activities  ( Con ' t ) 

Mr.  Jerome  Pereh 
Staff  Specialist  for  Materials 
and  Structures 
OUSDR&E ,  The  Pentagon 
loom  3D1089 

Washington,  D.C.  20301 

Chief,  Airframe  and  Equipment  Branch 
FS-120 

Office  of  Flight  Standards 
Federal  Aviation  Agency 
Washington,  D.C.  20553 

National  Academy  of  Sciences 
National  Research  Council 
Ship  Hull  Research  Committee 
2101  Constitution  Avenue 
Washington,  D.C.  20418 
Attn:  Mr.  A.  R.  Lytle 

National  Science  Foundation 
Engineering  Mechanics  Section 
Division  of  Engineering 
Washington,  D.C.  20550 

Picatinny  Arsenal 

Plastics  Technical  Evaluation  Center 
Attn:  Technical  Information  Section 
Dover,  New  Jersey  07801 

Maritime  Administration 
Office  of  Maritime  Technology 
14th  and  Constitution  Avenue,  NW. 
Washington,  D.C.  20230 

PART  2  -  Contractors  and  Other  Technical 
Collaborators 

Universities 

Dr.  J.  Tinsley  Oden 
University  of  Texas  at  Austin 
345  Engineering  Science  Building 
Austin,  Texas  78712 

Professor  Julius  Miklovitz 
California  Institute  of  Technology 
Division  of  Engineering 
and  Applied  Sciences 
Pasadena,  California  91109 


Universities  (Con't) 

Dr.  Harold  Liebovitz,  Dean 
School  of  Engineering  and 
Applied  Science 
George  Washington  University 
Washington,  D.C.  20052 

Professor  Eli  Sternberg 
California  Institute  of  Technology 
Division  of  Engineering  and 
Applied  Sciences 
Pasadena,  California  91109 

Professor  Paul  M.  Naghdi 
University  of  California 
Department  of  Mechanical  Engineering 
Berkeley,  California  94720 

Professor  A.  J.  Durelli 
Oakland  University 
School  of  Engineering 
Rochester,  Missouri  48063 

Professor  F.  L.  DiMaggio 
Columbia  University 
Department  of  Civil  Engineering 
New  York,  New  York  10027 

Professor  Norman  Jones 

The  University  of  Liverpool 

Department  of  Mechanical  Engineering 

P.  0.  Box  147 

Brown low  Hill 

Liverpool  L69  3BX 

England 

Professor  E.  J.  Skudrzyk 
Pennsylvania  State  University 
Applied  Research  Laboratory 
Department  of  Physics 
State  College,  Pennsylvania  16801 

Professor  J.  Rlosner 
Polytechnic  Institute  of  New  York 
Department  of  Mechanical  and 
Aerospace  Engineering 
333  Jay  Street 
Brooklyn,  New  York  11201 

Professor  R.  A.  Schapery 
Texas  ASM  University 
Department  of  Civil  Engineering 
College  Station,  Texas  77843 


Universities  (Con’t.) 
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Professor  Walter  D.  Pilkey 
Univeraity  of  Virginia 
Reaearch  Laboratories  for  the 
Engineering  Sciences  and 
Applied  Sciences 
Charlottesville,  Virginia  22901 

Profeaaor  K.  D.  Willmert 
Clarkson  College  of  Technology 
Department  of  Mechanical  Engineering 
Potsdam,  New  York  13676 

Dr.  Walter  E.  Haisler 
Texas  A&M  University 
Aerospace  Engineering  Department 
College  Station,  Texas  77843 

Dr.  Hussein  A.  (Camel 
University  of  Arizona 
Department  of  Aerospace  and 
Mechanical  Engineering 
Tucson,  Arizona  85721 

Dr.  S.  J.  Fenves 
Carnegie-Mellon  University 
Department  of  Civil  Engineering 
Schenley  Park 

Pittsburgh,  Pennsylvania  15213 

Dr.  Ronald  L.  Huston 
Department  of  Engineering  Analysis 
University  of  Cincinnati 
Cincinnati,  Ohio  45221 

Professor  G.  C.  M.  Sih 
Lehigh  University 
Institute  of  Fracture  and 
Solid  Mechanics 
Bethlehem,  Pennsylvania  18015 

Professor  Albert  S.  Kobayashi 
University  of  Washington 
Department  of  Mechanical  Engineering 
Seattle,  Washington  98105 

Professor  Daniel  Frederick 
Virginia  Polytechnic  Institute  and 
State  University 

Department  of  Engineering  Mechanics 
Blacksburg,  Virginia  24061 


Universities  (Con't) 

Professor  A.  C.  Eringen 
Princeton  University 
Department  of  Aerospace  and 
Mechanical  Sciences 
Princeton,  New  Jersey  08540 

Professor  E.  H.  Lee 

Stanford  University 

Division  of  Engineering  Mechanics 

Stanford,  California  94305 

Professor  Albert  I.  Ring 
Wayne  State  University 
Biomechanics  Research  Center 
Detroit,  Michigan  48202 

Dr.  V.  R.  Hodgson 
Wayne  State  University 
School  of  Medicine 
Detroit,  Michigan  48202 

Dean  B.  A.  Boley 
Northwestern  University 
Department  of  Civil  Engineering 
Evanston,  Illinois  60201 

Professor  P.  G.  Hodge,  Jr. 
University  of  Minnesota 
Department  of  Aerospace  Engineering 
and  Mechanics 

Minneapolis,  Minnesota  55455 

Dr.  D.  C.  Drucker 
University  of  Illinois 
Dean  of  Engineering 
Urbana,  Illinois  61801 

Professor  N.  M.  Nevmark 
University  of  Illinois 
Department  of  Civil  Engineering 
Urbana,  Illinois  61803 

Professor  E.  Reissner 
University  of  California,  San  Diego 
Department  of  Applied  Mechanics 
La  Jolla,  California  92037 

Professor  Willims  A.  Nash 
University  of  Massachusetts 
Department  of  Mechanics  and 
Aerospace  Engineering 
Amherst,  Massachusetts  01002 


Universities  (Con't) 
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Universities  (Con't) 


Professor  G.  Herrmann 
Stanford  University 
Departsient  of  Applied  Mechanics 
Stanford,  California  94305 

Professor  J.  D.  Achenbach 
Northwest  University 
Department  of  Civil  Engineering 
Evanston,  Illinois  60201 

Professor  S.  B.  Dong 
University  of  California 
Department  of  Mechanics 
Los  Angeles,  California  90024 

Professor  Burt  Paul 
University  of  Pennsylvania 
Tovne  School  of  Civil  and 
Mechanical  Engineering 
Philadelphia,  Pennsylvania  19104 

Professor  H.  W.  Liu 
Syracuse  University 
Department  of  Chemical  Engineering 
and  Metallurgy 
Syracuse,  New  York  13210 

Professor  S.  Bodner 
Technion  R6D  Foundation 
Haifa,  Israel 

Professor  Werner  Goldsmith 
University  of  California 
Department  of  Mechanical  Engineering 
Berkeley,  California  94720 

Professor  R.  S.  Rivlin 
Lehigh  University 
Center  for  the  Application 
of  Mathemat ics 

Bethlehem,  Pennsylvania  18015 

Professor  F.  A.  Cozzarelli 
State  University  of  New  York  at 
Buffalo 

Division  of  Interdisciplinary  Studies 
Karr  Parker  Engineering  Building 
Chemistry  Road 
Buffalo,’ New  York  14214 


Professor  Joseph  L.  Rose 
Drexel  University 

Department  of  Mechanical  Engineering 
and  Mechanics 

Philadelphia,  Pennsylvania  19104 

Professor  B.  K.  Donaldson 
University  of  Maryland 
Aerospace  Engineering  Department 
College  Park,  Maryland  20742 

Professor  Joseph  A.  Clark 
Catholic  University  of  America 
Department  of  Mechanical  Engineering 
Washington,  D.C.  20064 

Dr.  Samuel  B.  Batdorf 
University  of  California 
School  of  Engineering 
and  Applied  Science 
Los  Angeles,  California  90024 

Professor  Isaac  Fried 
Boston  University 
Department  of  Mathematics 
Boston,  Massachusetts  02215 

Professor  E.  Krempl 
Rensselaer  Polytechnic  Institute 
Division  of  Engineering 
Engineering  Mechanics 
Troy,  New  York  12181 

Dr.  Jack  R.  Vinson 
University  of  Delaware 
Department  of  Mechanical  and  Aerospa 
Engineering  and  the  Center  for 
Composite  Materials 
Newark,  Delaware  19711 

Dr.  J.  Duffy 
Brown  University 
Division  of  Engineering 
Providence,  Rhode  Island  02912 

Dr.  J.  L.  Swedlow 
Caraegie-Mel Ion  University 
Department  of  Mechanical  Engineering 
Pittsburgh,  Pennsylvania  15213 


Universities  (Con't) 
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Universities  (Con't) 


Dr.  V.  K.  Varadan 

Ohio  State  University  Research  Foundation 
Department  of  Engineering  Mechanics 
Columbus,  Ohio  43210 

Dr.  Z.  Hash  in 

University  of  Pennsylvania 
Department  of  Metallurgy  and 
Materials  Science 
College  of  Engineering  and 
Applied  Science 

Philadelphia,  Pennsylvania  19104 

Dr.  Jackson  C.  S.  Yang 
University  of  Maryland 
Department  of  Mechanical  Engineering 
College  Park,  Maryland  20742 

Professor  T.  Y.  Chang 
University  of  Akron 
Department  of  Civil  Engineering 
Akron,  Ohio  44325 

Professor  Charles  W.  Bert 
University  of  Oklahoma 
School  of  Aerospace,  Mechanical, 
and  Nuclear  Engineering 
Norman,  Oklahoma  73019 

Professor  Satya  N.  Atluri 
Georgia  Institute  of  Technology 
School  of  Engineering  and 
Mechanics 

Atlanta,  Georgia  30332 

Professor  Graham  F.  Carey 
University  of  Texas  at  Austin 
Department  of  Aerospace  Engineering 
and  Engineering  Mechanics 
Austin,  Texas  78712 

Dr .  S .  S .  Wang 
University  of  Illinois 
Deoartment  of  Theoretical  and 
Applied  Mechanics 
Urbana,  Illinois  61801 

Professor  J.  F.  Abel 
Cornell  University 
Department  of  Theoretical 
and  Applied  Mechanics 
Ithaca,  New  York  14853 


Professor  V.  H.  Neubert 
Pennsylvania  State  University 
Department  of  Engineering  Science 
and  Mechanics 

University  Park,  Pennsylvania  16802 

Professor  A.  W.  Leissa 
Ohio  State  University 
Department  of  Engineering  Mechanics 
Columbus,  Ohio  43212 

Professor  C.  A.  Brebbia 
University  of  California,  Irvine 
Department  of  Civil  Engineering 
School  of  Engineering 
Irvine,  California  92717 

Dr.  George  T.  Hahn 
Vanderbilt  University 
Mechanical  Engineering  and 
Materials  Science 
Nashville,  Tennessee  37235 

Dean  Richard  H.  Gallagher 
University  of  Arizona 
College  of  Engineering 
Tucson,  Arizona  85721 

Professor  E.  F.  Rybicki 
The  University  of  Tulsa 
Department  of  Mechanical  Engineering 
Tulsa,  Oklahoma  74104 

Dr.  R.  Haftka 

Illinois  Institute  of  Technology 
Department  of  Mechanics  and  Mechanics 
and  Aerospace  Engineering 
Chicago,  Illinois  60616 

Professor  J.  G.  de  Oliveira 
Massachusetts  Institute  of  Technology 
Department  of  Ocean  Engineering 
77  Massachusetts  Avenue 
Cambridge,  Massachusetts  02139 

Dr.  Bernard  W.  Shaffer 
Polytechnic  Institute  of  New  York 
Route  110 

Farmingdale,  New  York  11735 
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Industry  and  Research  Institutes 


Industry  and  Research  Institutes  (Con't) 


Dr.  Norman  Robbs 
Raman  AviDvne 
Division  of  Raman 
Sciences  Corporation 
Burlington,  Massachusetts  01803 

Argonne  National  Laboratory 
Library  Services  Department 
9700  South  Cass  Avenue 
Argonne,  Illinois  60440 

Dr.  M.  C.  Junger 
Cambridge  Acoustical  Associates 
54  Rindge  Avenue  Extension 
Cambridge,  Massachusetts  02140 

Mr.  J.  H.  Torrance 
General  Dynamics  Corporation 
Electric  Boat  Division 
Groton,  Connecticut  06340 

Dr.  J.  E.  Greenspon 

J.  G.  Engineering  Research  Associates 
3831  Menlo  Drive 
Baltimore,  Maryland  21215 


Dr.  T.  L.  Geers 

Lockheed  Missiles  and  Space  Company 

3251  Hanover  Street 

Palo  Alto,  California  94304 

Mr.  William  Cayvood 
Applied  Physics  Laboratory 
Johns  Hopkins  Road 
Laurel,  Maryland  20810 

Dr.  Robert  E.  Dunham 
Pacifica  Technology 
P.0.  Box  148 

Del  Mar,  California  92014 

Dr.  M.  F.  Ranninen 
Battelle  Columbus  Laboratories 
505  Ring  Avenue 
Columbus,  Ohio  43201 

Dr.  A.  A.  Hochrein 
Daedalean  Associates,  Inc. 
Springlake  Research  Road 
15110  Frederick  Road 
Woodbine,  Maryland  21797 


Newport  News  Shipbuilding  and 
Dry  Dock  Company 
Library 

Newport  News,  Virginia  23607 


Dr.  James  W.  Jones 
Swanson  Service  Corporat ion 
P.0.  Box  5415 

Huntington  Beach,  California  92646 


Dr.  W.  F.  Bozich 

McDonnell  Douglas  Corporation 

5301  Bolsa  Avenue 

Huntington  Beach,  California  92647 

Dr.  H.  N.  Abramson 
Southwest  Research  Institute 
8500  Culebra  Road 
San  Antonio,  Texas  78284 

Dr.  R.  C.  DeHart 

Southwes t  Research  Ins t it ut e 

8500  Culebra  Road 

San  Antonio,  Texas  78284 

Dr.  M.  L.  Baron 
Weidlinger  Associates 
110  East  59th  Street 
New  York,  New  York  10022 


Dr.  Robert  E.  Nickel  1 
Applied  Science  and  Technology 
3344  North  Torrey  Pines  Court 
Suite  220 

La  Jolla,  California  92037 

Dr.  Revin  Thomas 
Westinghouse  Electric  Corp. 
Advanced  Reactors  Division 
P.  0.  Box  158 

Madison,  Pennsylvania  15663 

Dr.  H.  D.  Hibbitt 
Hibbitt  &  Rarlsson,  Inc. 

132  George  M.  Cohan  Boulevard 
Providence,  Rhode  Island  02903 

Dr.  R.  D.  Mindlin 
89  Deer  Hill  Drive 
Ridgefield,  Connecticut  06877 
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Industry  and  Research  Institutes  (Con't) 

Dr.  Richard  E.  Dane 
Mega  Engineering 
11961  Tech  Road 

Silver  Spring,  Maryland  20904 

Mr.  G.  M.  Stanley 
Lockheed  Palo  Alto  Research 
Laboratory 
3251  Hanover  Street 
Palo  Alto,  California  94304 

Mr.  R.  L.  Cloud 

Robert  L.  Cloud  Associates,  Inc. 

2972  Adeline  Street 
Berkeley,  California  94703 


